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Discovery of a planet around the K giant star 4 UMa 

M.P. Dollinger 1 ' 5 , A. P. Hatzes 2 , L. Pasquini 1 , E. W. Guenther 2 , M. Hartmann 2 , L. Girardi 3 , and M. 

Esposito 2 ' 



o 
o 



(N 



> 

(N 

m 
o 

o 

a: 

6 

C3 



European Southern Observatory, Karl-Schwarzschild-Strasse 2, 85748, Garching bei Miinchen, Germany 

Thiiringer Landessternwarte Tautenburg, Sternwarte 5, D-07778 Tautenburg, Germany 

INAF-Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, F35122 Padova, Italy 

Dipartimento di Fisica "E.R. Caianiello", Universita di Salerno, via S. Allende, 84081 Baronissi (Salerno), Italy 

Max-Planck-Institut fur Astrophysik, Garching bei Miinchen, Germany 



Received; accepted 



ABSTRACT 



Context. For the past 3 years we have been monitoring a sample of 62 K giant stars using precise stellar radial velocity 
measurements with high accuracy taken at the Thiiringer Landessternwarte Tautenburg. 

Aims. To search for sub-stellar companions to giant stars and to understand the nature of the diverse radial velocity 
variations exhibited by K giant stars. 

Methods. We present precise stellar radial velocity measurements of the K1III giant star 4 UMa (HD 73108). These 
were obtained using the coude echelle spectrograph of 2-m Alfred Jensch Telescope. The wavelength reference for the 
radial velocity measurements was provided by an iodine absorption cell. 

Results. Our measurements reveal that the radial velocity of 4 UMa exhibits a periodic variation of 269.3 days with 
a semiamplitude K — 216.8 ms -1 . A Keplerian orbit with an eccentricity, e = 0.43 ± 0.02 is the most reasonable 
explanation for the radial velocity variations. The orbit yields a mass function, f(m) = (2.05 ± 0.24) x 10~ 7 M©. From 
our high resolution spectra we calculate a metallicity of —0.25 ± 0.05 and derive a stellar mass of 1.23 M© ± 0.15 for 
the host star. 

Conclusions. The K giant star 4 UMa hosts a substellar companion with minimum mass M sin i = 7.1 ± 1.6 Mj up iter- 
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1. Introduction 
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More than 200 extrasolar planets around main sequence 
stars have been detected via radial velocity method (RV), 
however most of these are around stars with masses of ap- 
proximately one solar mass. We thus have a very poor un- 
derstanding of how the stellar mass influences planet for- 
mation, particularly for more massive stars. Although a few 
radial velocity searches have been undertaken for planets 
around low mass stars (Delfosse et al. 1999; Endl et al. 
2003), there have been fewer attempts to search for plan- 
ets around intermediate mass stars. The reason is that RV 
searches are ill-suited for early-type, more massive main 
sequence stars. These stars are hotter and thus have fewer 
spectral lines, and these are broadened significantly by the 
high rotation rates found among early-type A-F stars. 

In spite of the lower RV precision, there are on-going 
attempts to search for planets around intermediate mass 
A-F stars (Galland et al. 2005a) and these have produced 
two substellar companions. A 9.1 Mj up i te r mass companion 
was found in a 388-d orbit around the F-type star (M = 
1.25 M dot) by Galland et al. (2005b). A candidate brown 
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dwarf (25 Mj up i te r) in a 28-day orbit was found around the 
A9V star HD 180777 (Galland et al. 2006). 

An alternative approach to search for planets around 
more massive stars is to look at intermediate mass stars 
that have evolved off the main sequence and up the giant 
branch. These stars are cooler and thus have a plethora 
of stellar lines. This along with their slower rotation rates 
make them amenable to high precision RV measurements. 
Two difficulties are encountered in this approach. First, 
unlike for main sequence stars, K giants of widely differ- 
ent masses can have similar effective temperatures. Second, 
giant stars show intrinsic variations due to stellar oscilla- 
tions. These short period (2-10 days) variations are likely 
caused by radial and/or nonradial p-mode oscillations (e.g. 
Hatzes & Cochran 1994). These oscillations add intrinsic 
RV "noise" making the detection of extrasolar planets more 
difficult. In spite of these challenges, searching for planets 
around giant stars can give us information about the pro- 
cess of planet formation around intermediate mass stars. 

The first indications of sub-stellar companions around 
giant stars was found by Hatzes & Cochran (1993) who dis- 
covered long period RV variations in three K giant stars. 
They proposed two viable hypotheses for these variations: 
sub-stellar companions or rotational modulation. The ex- 
pected rotational periods of K giants are several hundreds 
of days which are comparable to the observed RV periods. 
If a large surface inhomogeneity (e.g. starspot) were on the 
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surface this would create distortions of the spectral line pro- 
files which would be detected as an RV variation with the 
rotation period of the star. For these reasons the nature of 
the long period RV variations was not clear. 

Subsequent studies have established that giant stars can 
indeed host extrasolar planets. Frink et al. (2002) discov- 
ered the first extrasolar planet around the K giant star HD 
137759 (t Dra). This was followed by the discovery of sub- 
stellar companions to the K2 III star HD 47536 (Setiawan 
et al. 2003a) and HD 122430 (Setiawan 2003b). In the same 
year Sato et al. (2003) reported a planetary companion 
around G9 III HD 104985. Sub-stellar companions have 
also been reported for HD 11977 (Setiawan et al. 2005) 
and HD 13189 (Hatzes et al. 2005). More recently Hatzes 
et al. (2006) confirmed that the initial RV variations found 
by Hatzes & Cochran (1993) in /3 Gem were in fact due to 
a planetary companion. This was confirmed by Rcffcrt at 
al. (2006). 

For most of these discoveries the planet hypothesis for 
the RV variations was established due to a lack of variations 
in other measured quantities with the RV period. Spots or 
stellar pulsations are expected to also produce variations in 
cither activity indicators (Ca II H&K), photometry, and/or 
spectral line shape variations with the RV period. These 
variations were not found in the giant stars claimed to host 
extrasolar planets. In the case of l Dra the high eccentricity 
of the orbit established its Keplerian nature. 

2. Observations and data analysis 

4 UMa (= HD 73108 = HR 3403 = HIP 42527) be- 
longs to a sample of K giant stars that we have 
observed since February 2004 at the Thuringia State 
Observatory (Thuringer Landessternwarte Tautenburg or 
TLS). Observations were made using the 2m Alfred Jensch 
Telescope as part of the Tautenburg Planet Search Program 
(TOPS). The TOPS program uses the coude echelle spec- 
trograph which provides a resolving power, R = 67000 and 
a wavelength range of 4700-7400 Angstroms. A total of 
46 spectra (nightly averages) were taken of 4 UMa using 
the iodine cell. The total exposure time ranged between 

5 and 10 minutes depending on the weather conditions 
and this resulted in a typical signal-to-noise ratio typi- 
cally greater than 150. The strategy was to make obser- 
vations for this star 1-3 times per month, weather permit- 
ting. The standard CCD data reductions (bias-subtraction, 
flat-fielding and spectral extraction) were performed using 
reduction routines of the Image Reduction and Analysis 
Facility (IRAF). 

An iodine absorption cell placed in the optical path pro- 
vided the wavelength reference for the velocity measure- 
ments. The RVs were calculated by modeling the observed 
spectra with a high signal-to- noise ratio template of the star 
taken without the iodine cell and a scan of our iodine cell 
taken at very high resolution with the Fourier Transform 
spectrometer (FTS) of the McMath-Pierce telescope at Kitt 
Peak. The relative velocity shift between stellar and iodine 
absorption lines as well as the temporal and spatial vari- 
ations of the instrument profile were calculated for each 
observation. For more details about the spectrograph, RV 
program, and typical measurement precision see Hatzes et 
al. (2005). For the bright K giant stars in our program we 
typically achieve an RV accuracy of about 3-5 ms' 1 . A 
more detailed discussion of the RV accuracy for the K gi- 



ant stars will be presented in a forthcoming paper on the 
complete Tautenburg sample. 

3. The properties of the star 4 UMa 

4 UMa has a visual magnitude of V = 5.79 mag and is clas- 
sified in SIMBAD as K1III star. The Hipparcos parallax is 
12.92 ± 0.71 mas and this implies an absolute magnitude 
M v = 0.146 ± 0.119 mag. The stellar parameters of 4 UMa 
are summarized in Table 1. These were obtained either from 
the literature or were derived from our analysis of the stel- 
lar spectra taken without the iodine cell. Our high-quality 
spectra allowed us to determine accurate Fe abundances as 
well as the effective temperature, T c ff, the surface gravity, 
log g, and the microturbulcncc velocity, £. The results of 
this analysis for 4 UMa and the rest of the Tautenburg sam- 
ple will be presented in more detail in a forthcoming paper. 
The metallicity, T e fj, and the absolute V-band magnitude 
as derived from Hipparcos parallaxes were used as input 
values to estimate the mass, age, radius, (B-V)o and, in an 
alternative way, the surface gravity of the program stars by 
comparing these to theoretical isochrones and a modified 
version of J^rgensen & Lindegren's (2005) method. A de- 
tailed description about the procedure is given in da Silva et 
al. (2006). Previous investigations have also derived stellar 
parameters for 4 UMa. McWilliam (1990) obtained the fol- 
lowing values for the T c ff, surface gravity log g and [Fe/H]: 
4370 K, 2.45 and —0.26. The corresponding values by Luck 
(1991) were 4400 K, 1.61 and -0.20, respectively. Both of 
these values are in excellent agreement with our results. 

Table 1. Stellar parameters of 4 UMa 



Spectral type 


K1III 


HIPPARCOS 


mv 


5.787 ± 0.005 


mag 


M v 


0.146 ±0.119 


mag 


B-V 


1.197 ±0.005 


mag 


Parallax 


12.92 ±0.71 


mas 


Distance 


62.39 ± 3.43 


pc 


Mass (a) 


1.234 ±0.15 


M© 


p(<0 
-ft* 


18.11 ± 1.47 


R 


Age (a) 


4.604 ±2.0 


Gyr 


rp(a) 


4415 ± 70 


K 


[Fe/H] (a) 


-0.25 ± 0.04 


dex 


log<? (a) 


1.8 ±0.15 


dex 


micro turbulence < - a - > 


1.2 ±0.8 


kms -1 



Dollinger et al. in preparation 



4. Radial velocity variations and orbital solution 

The time series of our RV measurements for 4 UMa is 
shown in Figure 1 . There is an obvious sinusoidal variation 
in the RV curve with a period of approximately 270 days. 
This docs not appear to be a pure sine wave and is thus 
the first hint of Keplerian motion. 

An analysis using a Lomb-Scargle periodogram (Lomb 
1976, Scargle 1982) confirmed the presence of strong power 
at a frequency v = 0.0038 cd" 1 (P = 271 days). The false 
alarm probability (FAP) of this peak using the prescription 
in Scargle (1982) is estimated to be FAP w 10~ 8 . 
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Fig. 1. (Top) Radial velocity measurements for 4 UMa. The 
solid line is the orbital solution. (Bottom) Residual RV vari- 
ations after subtracting the contribution of the planet orbit. 

Table 2. Orbital parameters for the companion to 4 UMa. 



Period [days] 

T^periastron [JD] 

Kfms- 1 ] 

o-tO-CMms- 1 ] 

e 

w[deg] 

f(m) [solar masses] 
a[AU] 



269.3 ± 1.96 
52987.3936 ± 4.31 
215.55 ± 7.10 
28.8 

0.432 ± 0.024 

23.81 ± 4.42 

(2.045 ± 0.247) x 10~ 7 

0.87 ± 0.04 



An orbital solution using the RV data resulted in a pe- 
riod, P = 269.3 ± 1.96 days and an eccentricity e = 0.43 
± 0.02. The corresponding mass function is f(m) = (2.045 
± 0.247) x 10~ 7 M Q . Using our derived stellar mass of M* 
= 1.23 Mq ± 0.15 we calculated a minimum mass for the 
companion of M sin i = 7.1 ± 1.6 Mj up i te r- All the or- 
bital elements are listed in Table 2. Figure 2 shows the RV 
variations phase-folded to the orbital period. 

Although the orbital fit to the data is good, there are 
points that significantly from the solution. In fact, the rms 
scatter, c, about the orbital solution is rather large, a « 
30 ms' 1 , or about a factor of ten larger than our expected 
errors. We suspect that this results from another periodic 
signal in the data. The lower panel of Figure 1 shows the 
residual RV variations after removal of the orbital contri- 
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Fig. 2. Radial velocity measurements for 4 UMa phased to 
the orbital period. The line represents the orbital solution. 

bution due to the planetary companion. There arc clear 
variations with a period much larger than the orbital pe- 
riod. This is confirmed by the Lomb-Scargle periodogram of 
the radial velocity residuals shown in Figure 3. This shows 
a strong peak at a frequency of v = 0.00156 cd _1 (P = 
641 days). The statistical significance of this signal was es- 
timated using a 'bootstrap randomization' technique. The 
measured RV values were randomly shuffled keeping the 
observed times fixed and a periodogram for the shuffled 
data computed. After 2x 10 5 shuffles there was no instance 
where the random fake data periodogram data had higher 
power than the data periodogram in the frequency range. 
This indicates that the FAP < 2 x 10~ 6 . This signal is thus 
statistically significant. 

We also note that K giants show RV variations on time 
scales of several days and amplitudes of up to 100 ms _1 . 
This may also contribute to the scatter of the RV points 
about the orbital solution. 



5. Hipparcos photometry 

We also analyzed the Hipparcos photometry of 4 UMa in or- 
der to see if variations with the 270 day period are present. 
Hipparcos obtained 27 observations ('daily" averages') over 
a 3 year period. Figure 4 shows the Lomb-Scargle peri- 
odogram of the Hipparcos photometry after averaging the 
multiple measurements taken within a few hours of each 
other. There was on obvious outlier that had a value more 
than lOtr from the mean that was eliminated from the data 
prior to calculating the periodogram. Although there is a 
weak peak near the orbital frequency this is not significant 
and there are several peaks with more power. The highest 
peak is at a frequency of 0.0252 cd _1 (P = 39.698 days). 
Figure 5 shows the photometry phased to this period. The 
false alarm probability was estimated using the bootstrap 
randomization technique. After 100,000 such shuffles 19% of 
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Fig. 3. Lomb-Scargle periodogram of the RV residuals of 4 
UMa. The strong peak corresponds to a period of 641 days. 



6. Discussion 

Our radial velocity measurements indicate that the giant 
star 4 UMa hosts a giant extrasolar planet. We are confi- 
dent that these RV variations are due to a sub-stellar com- 
panion and not rotational modulation for three reasons. 
First, the orbit is highly eccentric and such a saw-toothed 
RV pattern is difficult to reproduce with rotational modu- 
lation from spots or possibly even from stellar oscillations. 
Indeed, the high eccentricity of the companion to i Dra 
was the convincing argument that this was a true compan- 
ion. Second, the Hipparcos photometry does not show any 
significant variations at the orbital period. Finally, our pe- 
riodogram analysis shows a statistically significant period 
at 641 days in the RV residuals for 4 UMa. We have de- 
tected two periods of several hundred days; both cannot 
be due to rotation. Because the RV variations of the dom- 
inant 269-d RV period is well-fit by an eccentric Kcplerian 
orbit, we believe that this signal is due to a planetary com- 
panion. The 641-day period we have detected may indeed 
be due to rotational modulation, although a similar signal 
does not appear in the Hipparcos photometry. We caution 
the reader, however, that the Hipparcos photometry was 
not contemporaneous with our RV measurements. 



the random data periodograms showed power greater than 
the real periodogram. This peak is most likely not signifi- 
cant. It it were real it be could be due to stellar oscillations. 
For example, Hatzes & Cochran (1998) detected a similar 
period (« 50 days) in the spectral line bisector variations 
in Aldebaran which they attributed to oscillations. 



s 

o 





.04 



.05 



.02 .03 
Frequency (cycles/day) 

Fig. 4. Lomb-Scargle periodogram of the Hipparcos pho- 
tometry for 4 UMa. The horizontal line marks the orbital 
frequency. 
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Fig. 5. Hipparcos photometry for 4 UMa phased to the 39.6 
day photometric period. 

At this point we do not wish to speculate as to the 
nature of the 641-day period that we have detected in the 
RV residuals. It could well be due to rotational modulation, 
or possibly due to a second companion. We are continuing 
to monitor this star and are currently making a detailed 
analysis of other observed quantities (Ha for example) to 
see if this in fact the rotation period of the star. 

The high eccentricity of the orbit with a value of e = 
0.43 ± 0.02 is not unusual for planet hosting stars, including 
giant stars. The companion to t Dra has the highest eccen- 
tricity (e = 0.7) for a planet around a giant star. Relatively 
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high eccentricities have also been found for the companions 
to HD 11977 (e = 0.4, Setiawan et al. 2005) and HD 13189 
(e = 0.27, Hatzes et al. 2006). It seems that the planetary 
companions to giant stars can have the wide range of or- 
bital eccentricities that are shown by planets orbiting main 
sequence stars. 

The abundance analysis for 4 UMa shows that this star 
is slightly metal poor ([Fe/H] = -0.25 ± 0.05). This in- 
teresting because this the main sequence host stars of exo- 
planets tend to to be metal rich compared to stars that do 
not posses exoplanets (Santos et al. 2004). Other authors 
(Schuler et al. 2006; da Silva et al. 2006) have also found 
evidence that planet hosting giant stars are metal poor. 
However, given the small number of planets around giant 
stars it is premature to conclude that these contradict the 
planet-metallicity trend for main sequence stars. An anal- 
ysis of a larger sample of giant stars may show that those 
have higher metallicity indeed have a higher frequency of 
sub-stellar companions. 

The number of exoplanets around evolved stars is so far 
very limited. The discovery of the companion to 4 UMa 
adds to the growing list of planets evolved stars hosting 
sub-stellar companions (see Hatzes et al. 2006 and refer- 
ences therein). The estimated masses for the giant stars 
hosting planets range from 1 to 3 solar masses. In terms of 
the properties of planets around giant stars they have sim- 
ilar characteristics, masses in the range of of 2-14 Mj up i t er 
and orbital periods of several hundred days (see discus- 
sion by Hatzes et al. 2006). In other words, similar to the 
properties of the planet around 4 UMa. These discoveries 
are important because they probe a different stellar mass 
regime different from the main sequence objects that are 
the targets of most planet searches. 
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